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Abstract 

This paper investigates Starobinsky's model of inflation driven by the trace anomaly 
of conformally coupled matter fields. This model does not suffer from the problem of 
contrived initial conditions that occurs in most models of inflation driven by a scalar field. 
The universe can be nucleated semi-classically by a cosmological instanton that is much 
larger than the Planck scale provided there are sufficiently many matter fields. There are 
two cosmological instantons: the four sphere and a new "double bubble" solution. This 
paper considers a universe nucleated by the four sphere. The AdS/CFT correspondence is 
used to calculate the correlation function for scalar and tensor metric perturbations during 
the ensuing de Sitter phase. The analytic structure of the scalar and tensor propagators is 
discussed in detail. Observational constraints on the model are discussed. Quantum loops 
of matter fields are shown to strongly suppress short scale metric perturbations, which 
implies that short distance modifications of gravity would probably not be observable in the 
cosmic microwave background. This is probably true for any model of inflation provided 
there are sufficiently many matter fields. This point is illustrated by a comparison of 
anomaly driven inflation in four dimensions and in a Randall- Sundrum brane- world model. 
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1 Introduction 



Inflation |]J in the very early universe seems the only natural explanation of many observed 
features of our universe, in particular the recent measurements of a Doppler peak in the cosmic 
microwave background fluctuations ||. However, while it provides an appealing explanation for 
several cosmological problems, it provokes the natural question of why the conditions were such 
as to start inflation in the first place. 

The new inflationary scenario || [| was proposed primarily to overcome the problem of 
obtaining a natural exit from the inflationary era. In this model, the value of the scalar is 
supposed to be initially confined to zero by thermal effects. As the universe expands and cools 
these effects disappear, leaving the scalar field miraculously exposed on a mountain peak of the 
potential. If the low temperature potential is sufficiently flat near = then slow roll inflation 
will occur, ending when the field reaches its true minimum <p c . This scenario seems implausible 
because a high temperature would confine only the average or expectation value of the scalar to 
zero. Rather than be supercooled to a state with ~ locally, the field fluctuates and rapidly 
forms domains with cf) near ±0 C The dynamics of the phase transition is governed by the growth 
and coalescence of these domains and not by a classical roll down of the spatially averaged field 
<p ||. Because this and other problems, new inflation was largely abandoned in favor of chaotic 
inflation || in which it is just assumed that that the scalar field was initially displaced from 
the minimum of the potential. One attempt to explain these initial conditions for inflation in 
terms of quantum fluctuations of the scalar field seems to lead to eternal inflation at the Planck 
scale 0, at which the theory breaks down. Another attempt, using the Hartle-Hawking "no 
boundary" proposal ||, found that the most probable universes did not have enough inflation 
||. No satisfactory answer to the question of why the scalar field was initially displaced from 
the minimum of its potential has been found. 

In this paper we will reconsider an earlier model, in which inflation is driven by the trace 
anomaly of a large number of matter fields. The Standard Model of particle physics contains 
nearly a hundred fields. This is at least doubled if the Standard Model is embedded in a 
super symmetric theory. Therefore there were certainly a large number of matter fields present 
in the early universe, so the large N approximation should hold in cosmology, even at the 
beginning of the universe. In the large N approximation, one performs the path integral over 
the matter fields in a given background to obtain an effective action that is a functional of the 
background metric: 

exp(-W[g}) = J d[<f>] exp(-S[0; g]). (1.1) 

One then argues that the effect of gravitational fluctuations is small in comparison to the large 
number of matter fluctuations. Thus one can neglect graviton loops, and look for a stationary 
point of the combined gravitational action and the effective action for the matter fields. This is 
equivalent to solving the Einstein equations with the source being the expectation value of the 
matter energy momentum tensor: 

R ij -^Rg ij = torG{T ij ), (1.2) 

where 

<T«> = (1.3) 
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Finally, one can calculate linearized metric fluctuations about this stationary point metric and 
check they are small. This is confirmed observationally by measurements of the cosmic mi- 
crowave background, which indicate that the primordial metric fluctuations were of the order 



of 10" 5 



Matter fields might be expected to become effectively conformally invariant if their masses 
are negligible compared to the spacetime curvature. Classical conformal invariance is broken 



at the quantum level [JD]] (see [12|, [13| for reviews), leading to an anomalous trace for the 
energy-momentum tensor: 

//'•'(v;,> / o. (i.4) 

This trace is entirely geometrical in origin and therefore independent of the quantum state. In 
a maximally symmetric spacetime, the symmetry of the vacuum implies that the expectation 
value of the energy momentum tensor can be expressed in terms of its trace 

(0\Tu\0) = ^ w <o|r w |o>. (1.5) 

Thus the trace anomaly acts just like a cosmological constant for these spacetimes. Hence a 
positive trace anomaly permits a de Sitter solution to the Einstein equations H . 



This is very interesting from the point of view of cosmology, as pointed out by Starobinsky 
P~5| ] . Starobinsky showed that the de Sitter solution is unstable, but could be long-lived, and 
decays into a matter dominated Friedman-Robertson- Walker (FRW) universe. The purpose of 
Starobinsky's work was to demonstrate that quantum effects of matter fields might resolve the 
Big Bang singularity^. From a modern perspective, it is more interesting that the conformal 
anomaly might have been the source of a finite but significant period of inflation in the early 
universe. This inflation would be followed by particle production and (p)reheating during the 
subsequent matter dominated phase. Starobinsky's work is reviewed and extended by Vilenkin 
in jr7|. For a more recent discussion of the Starobinsky model, see [IS . 



Starobinsky showed that the de Sitter phase is unstable both to the future and to the past, so 
it was not clear how the universe could have entered the de Sitter phase. However, this problem 
can be overcome by an appeal to quantum cosmology, which predicts that the de Sitter phase 
of the universe is created by semi-classical tunneling from nothing. This process is mediated 



by a four sphere cosmological instanton [|H]]. One of the results of this paper is that the four 
sphere is not the only cosmological instanton in this model. 

In order to test the Starobinsky model, it is necessary to compare its predictions for the 
fluctuations in the cosmic microwave background (CMB) with observation. This was partly 
addressed by Vilenkin [I7|. Using an equation derived by Starobinsky |]I5|| , Vilenkin showed 



that the amplitude of long wavelength gravitational waves could be brought within observational 
limits at the expense of some fine-tuning of the coefficients parameterizing the trace anomaly. 
Density perturbations were discussed by Starobinsky in [^ . 

The analysis of Starobinsky and Vilenkin was complicated by the fact that tensor pertur- 
bations destroy the conformal flatness of a FRW background, making the effective action for 
matter fields hard to calculate. However, we now have a way of calculating the effective action 
for a particular theory, namely M = 4 U(N) super Yang-Mills theory, using the AdS/CFT 

1 Another paper Jl6| which discussed the effects of the trace anomaly in cosmology failed to obtain non-singular 
solutions because it included a contribution from a classical fluid. 
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correspondence pl| . In this paper we will calculate the effective action for this theory in a 
perturbed de Sitter background. This enables us to calculate the correlation function for metric 
perturbations around the de Sitter background. We can then compare our results with observa- 
tions. The fact that we are using the M = 4 Yang-Mills theory is probably not significant, and 
we expect our results to be valid for any theory that is approximately massless during the de 
Sitter phase. One might think that our results could shed light on the effects of matter interac- 
tions during inflation since AdS/CFT involves a strongly interacting field theory. However, as 
we shall explain, our results are actually independent of the Yang-Mills coupling. 

Our calculations will be performed in Euclidean signature (on the four sphere), and then 
analytically continued to Lorentzian de Sitter space. The condition that all perturbations are 
regular on the four sphere defines the initial quantum state for Lorentzian perturbations. The 
four sphere instanton is much larger than the Planck scale (since we are dealing with a large N 
theory), so there is a clear cut separation into background metric and fluctuation. 

We shall include in our action higher derivative counterterms, which arise naturally in the 
renormalization of the Yang-Mills theory. There are three independent terms that are quadratic 
in the curvature tensors: the Euler density, the square of the Ricci scalar and the square of the 
Weyl tensor. The former just contributes a multiple of the Euler number to the action. Metric 
perturbations do not change the Euler number, so this term has no effect. The square of the 
Ricci scalar has the important effect of adjusting the coefficient of the V 2 i? term in the trace 
anomaly. It is precisely this term that is responsible for the Starobinsky instability, so by varying 
the coefficient of the R 2 counter term we can adjust the duration of inflation. The Weyl-squared 
counterterm does not affect the trace anomaly but it can contribute to suppression of tensor 
perturbations. The effects of this term were neglected by Starobinsky and Vilenkin. They also 
neglected the effects of the non-local part of the matter effective action. We shall take full 
account of all these effects. 

Vilenkin showed that the initial de Sitter phase is followed by a phase of slow-roll inflation 
before inflation ends and the matter-dominated phase begins. Since the horizon size grows 
significantly during this slow-roll phase, it is important to investigate whether modes we observe 
today left the horizon during the de Sitter phase or during the slow-roll phase. If the present 
horizon size left during the de Sitter phase, we find that the amplitude of metric fluctuations 
can be brought within observational bounds if N, the number of colours, is of order 10 5 . Such 
a large value for N is rather worrying, which leads us to the second possibility, that the present 
horizon size left during the slow-roll phase. Our results then suggest that the coefficient of 
the R 2 term must be at most of order 10 8 , and maybe much lower, but N is unconstrained 
(except by the requirement that the large N approximation is valid so that AdS/CFT can be 
used). We also find that the tensor perturbations can be suppressed independently of the scalar 
perturbations by adjusting the coefficient of the Weyl-squared counterterm in the action. 

Inflation blows up small scale physics to macroscopic scales. This suggests that inflation 
may lead to observational consequences of small-scale modifications of Einstein gravity, such as 
extra dimensions. However, we find that the non-local part of the matter effective action has the 
effect of strongly suppressing tensor fluctuations on very small scales, a result first noted in flat 
space by Tomboulis [52]. This suggests that any small-scale modifications to four dimensional 
Einstein gravity would be unobservable in the CMB since matter fields would dominate the 
graviton propagator at the scales at which such modifications might be expected to become 
important. This result is probably not restricted to trace anomaly driven inflation since it is 
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simply a consequence of the presence of a large number of matter fields. As we have mentioned, 
there really are a large number of matter fields in the universe and these will suppress small-scale 
graviton fluctuations in any model of inflation. 

We illustrate this point by considering a Randall-Sundrum (RS) |£J version of the Starobin- 
sky model. In the RS model, our universe is regarded as a thin domain wall in anti-de Sitter 
space (AdS). RS showed that linearized four dimensional gravity is recovered on the domain wall 
at distances much larger that the AdS radius of curvature, but gravity looks five dimensional 
at smaller scales. Therefore, if the AdS length scale is taken to be small, then the RS model 
is a short distance modification of four dimensional Einstein gravity. We shall show that when 
the large N field theory is included, the effects of the matter fields dominate the RS correc- 
tions to the graviton propagator and render them unobservable. This work is an extension of 
our previous paper |24| to include the effects of scalar perturbations and the higher derivative 
counterterms in the action. 

This paper is organized as follows. We start in section by showing that the Starobinsky 
model has two instantons: the round four sphere and a new "double bubble" instanton. We 
consider only the four sphere instanton in this paper. In section [| we use the AdS / CFT cor- 
respondence to calculate the effective action of the large iV Yang-Mills theory on a perturbed 
four sphere. Coupling this to the gravitational action then allows us to compute the scalar and 
tensor graviton propagators on the four sphere. In section |], we discuss the analytic structure 
of our propagators. The tensor propagator is shown to be free of ghosts. In section |5|, we show 
how our Euclidean propagators are analytically continued to Lorentzian signature. Section |] 
discusses two observational constraints on the Starobinsky model, namely the duration of infla- 
tion and the amplitude of perturbations. In section [7], we use the RS version of the Starobinsky 
model as an example to illustrate how matter fields strongly suppress metric perturbations on 
small scales. Finally, we summarize our conclusions and suggest possible directions for future 
work. 



2 0(4) Instantons 
2.1 Introduction 

Homogeneous isotropic FRW universes are obtained by analytic continuation of cosmological in- 
stantons invariant under the action of an 0(4) isometry group. In other words, we are interested 
in instantons with metrics of the form 

ds 2 = da 2 + b{a) 2 dtt 2 3 . (2.1) 

We shall restrict ourselves to instantons with spherical topology, for which fe(er) vanishes at a 
"North pole" and a "South pole". Regularity requires that b'(a) = ±1 at these poles. (Instan- 
tons with topology R x S 3 also exist but are excluded by the no boundary proposal because 
they are non-compact.) The scale factor b(a) is determined by Einstein's equation^ 

Gij = 87rG/(Ty>, (2.2) 

2 We use a positive signature metric and a curvature convention for which a sphere or de Sitter space has 
positive Ricci scalar. 
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where the right hand side involves the expectation value of the energy momentum tensor of the 
matter fields, which we are assuming to come from the jV = 4 U (N) super Yang-Mills theory. 
(Tij) can be obtained for the most general quantum state of the Yang-Mills theory consistent 
with 0(4) symmetry by using the trace anomaly and energy conservation, as we shall describe 
below. 



2.2 The trace anomaly 

The general expression for the trace anomaly of our strongly coupled large N CFTQ was calcu- 
lated using AdS/CFT in [25]. It turns out that it is exactly the same as the one loop result for 
the free theory, which is given for a general CFT by the following equation [|1^, |H| 

ft fa) = cF-aG + dV 2 R (2.3) 

where F is the square of the Weyl tensor: 

F = C ijkl C^ kl , (2.4) 

G is proportional to the Euler density: 

G = R ijk iR ijkl - 4RijR ij + R 2 , (2.5) 

and the constants a, c and d are given in terms of the field content of the CFT by 

1 



360(4tt) 2 
1 

120(4vr) 2 



(N s + UN F + Q2N V ) , (2.6) 
(N s + 6N F + 12N V ) , (2.7) 



d= j^- 2 (N S + W F - UNy) , (2.8) 

where Ns is the number of real scalar fields, N F the number of Dirac fermions and Ny the 
number of vector fields. The coefficients a and c are independent of renormalization scheme but 
d is not. We have quoted the result given by zeta- function regularization or point-splitting; the 
result given by dimensional regularization has +12 instead of —18 as the coefficient of Ny JT2[ . 
In fact, d can be adjusted to any desired value by adding the finite counter term 

Set = ^ / d\^-gR 2 . (2.9) 

This counter term explicitly breaks conformal invariance. a is a dimensionless constant. The 
field content of the Yang-Mills theory is N s = QN 2 , N F = 2N 2 (there are AN 2 Majorana 
fermions, which is equivalent to 2N 2 Dirac fermions) and Ny = N 2 . This gives 

N 2 

a = c=—, d = 0. (2.10) 



3 We shall often refer to the Af = 4 Yang-Mills theory as a CFT even though it is not conformally invariant 
on the four sphere. 
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We have used the coefficient —18 for Ny when calculating d - this is the value predicted by 
AdS/CFT |[25|| . If d = then inflation never ends in Starobinsky's model. We shall therefore 
include the finite counter term, which does not change a or c but gives 

aiV 2 

rf=T7r^. (2.11) 
lo^ 

When we couple the Yang-Mills theory to gravity, the presence of Set implies that we are 
effectively dealing with a higher derivative theory of gravity. It is, of course, arbitrary whether 
one regards S ct as part of the gravitational action or as part of the matter action. We have 
adopted the latter perspective and therefore included an explicit factor of N 2 in the action 
(since there are 0(N 2 ) fields in the Yang-Mills theory). 



2.3 Energy conservation 



Having obtained the trace of the energy-momentum tensor, we can use energy-momentum con- 
servation to obtain the full energy-momentum tensor. Introduce the energy density p and 
pressure p, defined in an orthonormal frame by 

(T aa ) = -p, (T a p)=p5 aP . (2.12) 

The minus sign in the first expression arises because we are considering Euclidean signature. 
These must obey 

-p + 3p=(T), (2.13) 
and we also have the energy-momentum conservation equation 



P' + yKP + P) = 0. 

Eliminating p gives an equation for p: 

(b A p)' = -b 3 b'(T). 
Substituting in the expression for (T) and integrating gives 
3iV 2 



P 



8tt 2 6 4 



4 V 2 2 



(2.14) 



(2.15) 



(2.16) 



The expression for p is easily determined from equation 2.13 . The appearance of the constant 
of integration C shows that the quantum state can contain an arbitrary amount of radiation. 
Setting C = a/2 reproduces the energy-momentum tensor for the vacuum state. The cosmology 
resulting from the trace anomaly in the presence of an arbitrary amount of null radiation was 

The cosmological solutions obtained were generically singular. However, 
showed if this null radiation is not present (i.e., if C = a/2) then non-singular 



investigated in fl6 
Starobinsky 
solutions can be obtained. 

To conclude, we have found the energy-momentum tensor for a strongly coupled large iV 
Yang-Mills theory in the most general quantum state that is consistent with 0(4) symmetry. 
The effects of strong coupling do not show up in our energy-momentum tensor, which is of the 
same form as used in iTBI, ITBl. In the next subsection we shall use this result in the Einstein 



equations to determine the shape of the instanton. 
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2.4 Shape of the instanton 

Taking the aa component of the Einstein equation gives 



k/2 



Gr 



Substituting in our result for p gives 



b 2 



N 2 G 



b 2 



71 



1 



a 



'b'b' r 



-8nGp. 



b" 2 b' 2 b" 
262 + 



3b' 4 b' 2 



+ 



c 
b 1 



(2.17) 



(2.18) 



Regularity at the poles of the instanton requires b' — > ±1 as b — > 0. Substituting this into 
equation |2.18| , one finds that b" = and C = a/2 are also required for regularity at the poles. 
In other words, the no boundary proposal has singled out a particular class of quantum states 
for us, namely those that do not contain any radiation. These are precisely the states that can 
give rise to non-singular cosmological solutions. In our picture this is because such cosmological 
solutions can be obtained from a Euclidean instanton. 

It is convenient to introduce a length scale R defined by 



R 2 

We can now define dimensionless variables 

a = cr/R, 



N 2 G 

47T 



f(a) = b(a)/R. 



Equation [2.18| becomes 

1 - f' 2 (1 - f 2 ) 2 



P 



+ 2a 



2fT 
f 2 



f" 2 J' 2 f" 
f 2 P 



'f 



f 



12 



3 \h) +2^T + -T7 



./ 



P P 



(2.19) 



(2.20) 



(2.21) 



The boundary conditions at the poles are / = 0, /' = ±1, /" = (where a prime now denotes 
a derivate with respect to a). One solution to equation [2.21| is 



sin cr, 



(2.22) 



which simply gives us a round four sphere instanton. Note that the expression multiplying a 
vanishes for this solution. Another simple solution is 



(2.23) 



i.e. flat Euclidean space. 

In order to integrate p.21| numerically, we assume that a = is a regular "North pole" of 
the instanton. We start the integration at a — e. The boundary conditions for the integration 
are 

' (2.24) 



/( e ) = e + -f(0) e 3 + ... 
f'(e) = 1 + \p" (0)e 2 + . . . 



(2.25) 
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12 3 4 

Figure 1: Scale factor f(a) for a regular "double bubble" instanton with a = — 1 and f"'(0) = 
-2.05. 

f(e) = f'(0) e + ... (2.26) 

We shall neglect the higher order terms (denotes by the ellipses) in our numerical integration. 
It is important to retain all of the terms displayed in order to obtain f'"(e) = /"'(0) + . . . from 
the equation of motion. Note that f"'(0) is a free parameter. Our strategy is to choose the 
value of /"'(0) so that the instanton is compact and closes off smoothly at the South pole. 

The instanton is non-compact when f"'(0) > 0. The solution is flat Euclidean space when 
/'"(0) = 0. We shall therefore concentrate on f"'(0) < 0. The four sphere solution has f"'(0) = 
— 1. It is convenient to discuss the cases a > and a < separately. 

If a > then there are two types of behaviour, (i) —1 < /"'(0) < the instanton is non- 
compact. For f"'(0) close to —1, the scale factor increases to a local maximum and then starts 
to decrease. However, before reaching / = 0, the scale factor turns around again and increases 
indefinitely, (ii) f" r (0) < — 1. These instantons are compact but do not have a regular South 
pole since b' diverges there. They are the analogues of the singular instantons discussed in ||. 

If a < then there are two types of behaviour, (i) —1 < /"'(0) < 0. These instantons are 
compact with an irregular South pole, (ii) /"'(0) < —1. The scale factor of these instantons 
increases to a local maximum, decreases to a local minimum, then has another maximum before 
decreasing to zero at the South pole, which is irregular. The instanton therefore has two "peaks" . 
There is a critical value 7(a) such that for 7 < f"'(0) < — 1 the larger peak is near the North 
pole while for f"'(0) < 7, the larger peak is near the South pole. It follows that when f"'(0) = 7 
the peaks have the same size and the instanton is symmetrical about its equator with a regular 
South pole. The scale factor is shown in figure [l|. 

To summarize, if a < then there are two regular compact instantons, namely the round 
four sphere and a new "double bubble" instanton. We shall not have much to say about the new 
instanton in this paper since the lack of an analytical solution makes dealing with perturbations 
of this instanton rather difficult. 
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2.5 Analytic continuation 



The four sphere instanton can be analytically continued to Lorentzian signature by slicing at 
the equator a = n/2 and writing 

a = - - it/R, (2.27) 

which yields the metric on a closed de Sitter universe: 

ds 2 = -dt 2 + R 2 cosh 2 (t/R)dQl (2.28) 

The Hubble parameter is R~ l , which is much smaller than the Planck mass because N is large. 
A change of coordinate takes one from a closed FRW metric to an open FRW metric. 

The double bubble instanton can be analytically continued across its "equator" to give a 
closed FRW universe. Numerical studies suggest that this universe rapidly collapses. How- 
ever, this instanton can also be continued to an inflationary open universe (the details of the 
continuation are the same as in ||) and therefore may give rise to realistic cosmology. 



3 Metric perturbations 

3.1 Scalars, vectors and tensors 

In this section we shall calculate correlation functions for metric perturbations around our four 
sphere instanton. These can then be analytically continued to yield correlation functions in de 
Sitter space. The metric on the perturbed four sphere can be written 

ds 2 = (R 2A fij + hij)dx l dx 3 , (3.1) 

where 7^ denotes the metric on a unit four sphere. The perturbation can be decomposed into 
scalar, vector and tensor parts with respect to the four sphere: 

hij(x) = %(x) + 2V (i Xj)(x) + ViV j( p(x) + %i/;(x). (3.2) 

The connection on the unit four sphere is denoted V. % is a transverse traceless symmetric 
tensor with respect to the four sphere: 

V^ ij = 0\ = 0, (3.3) 
where indices i, j are raised and lowered with 7^. \% is a transverse vector: 

ViX* = 0. (3.4) 

There is a small ambiguity in our decomposition - it is invariant under 0— > <p + Y, ip — > ip + XY 
where Y satisfies 

VtVjY + \%Y = 0. (3.5) 

This equation can only be solved when A = 1. The solutions are simply the regular p = 1 
spherical harmonics on S* 4 , i.e., the regular p = 1 solutions of 

(V 2 +p(p + 3))Y = 0. (3.6) 
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The spherical harmonics are labelled with integers p,k,l,m with < \m\ < I < k < p. Hence 
there are five independent spherical harmonics withp = 1, given in terms of spherical harmonics 
Y\dm on the three sphere by 

smpYum, cos pY 0Q0 (3.7) 

where p is the polar angle on the four sphere. These five harmonics correspond to gauge 
transformations involving the five conformal Killing vector fields on the four sphere [p6|| . If we 



assume that if) is regular on S 4 then we can expand it in terms of spherical harmonics. We shall 
fix the residual gauge ambiguity by demanding that if} contain no contribution from the p = 1 
harmonics. 

It is possible to gauge away and x 1 through a coordinate transformation on the four sphere 
of the form x l — > x % — rf — d l rj, where rf is a transverse vector and r\ is a scalar. For the moment 
we shall use a general gauge but later we will assume that and \ % vanish. 

3.2 Matter effective action 

We need to calculate the action for metric perturbations. The hardest part to calculate is 
the effective action for the matter fields. This can be expanded around a round four sphere 
background: 

W = W^~^Jd 4 x^/(T l:j (x))h^(x) 

+ ~J d 4 x^ J d A x'^h i \x)(T ij (x)T kl (x'))h kl (x') + ... (3.8) 

Here 7 denotes the determinant of the metric on the sphere. If we know the one and two 
point function of the CFT energy momentum tensor on a round S 4 then we can calculate the 
effective action to second order in the metric perturbation. The one point function is given 
by the conformal anomaly on the round four sphere. In flat space, the 2-point function is 
determined entirely by conformal invariance. On the sphere, symmetry determines the 2-point 
function only up to a single unknown function [pTfl . However, the sphere is conformally flat so 
one can calculate the 2-point function on the sphere using a conformal transformation from flat 
space. The energy-momentum tensor transforms anomalously, so there will be a contribution 
from the trace anomaly in the transformation. Therefore, the 2-point function on the sphere is 
determined by two quantities, namely the 2-point function in flat space, and the trace anomaly. 
For the M = 4 super Yang-Mills theory that we are considering, both of these quantities are 
independent of the Yang-Mills coupling. It follows that the 2-point function on the sphere 
(or any other conformally flat space) must be independent of coupling. Therefore the effective 
action will be independent of coupling to second order in the metric perturbation so the effects 
of strong coupling will not show up in our results. 

For the moment, we shall consider the four sphere to have arbitrary radius R rather than 
using the value given by equation |2.19| . Introduce a fictional ball of AdS that has the sphere as 



its boundary. Let I, G be the AdS radius and Newton constant of this region. If we take I to 
zero then the sphere is effectively at infinity in AdS so we can use AdS/CFT to calculate the 
generating functional of the CFT on the sphere. In other words, I is acting like a cut-off in the 
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CFT and taking it to zero corresponds to removing the cut-off. However the relation 



P 2N 2 , . 

Cr 7T 



implies that if I is taken to zero then we must also take G to zero since N is fixed (and large). 

The CFT generating functional is given by evaluating the action of the bulk metric g that 
matches onto the metric h of the boundary j3(| and adding surface counterterms to cancel 
divergences as I, G — > 



W[h] = S EH [g] + S GH [g] + Si[h] + S 2 [h] + S 3 [h] + S ct [h], (3.10) 

where S^h denotes the five dimensional Einstein-Hilbert action with a negative cosmological 
constant: 

Seh ~- 



(3 - n) 



167 

the overall minus sign arises because we are considering a Euclidean signature theory. The 
second term in the action is the Gibbons-Hawking boundary term [B{J: 

Sgh = ~~ I d^xVhK, (3.12) 

07TCx J 

where K is the trace of the extrinsic curvature of the boundary and h the determinant of the 
induced metric. The first two surface counterterms are 

S ' = 8^/^ (313) 
S 2 = — 1 -~ [ d 4 xVhR, (3.14) 



32ttG 

where R now refers to the Ricci scalar of the boundary metric. The third counterterm isQ 

S 3 = (log(l/R) -P)J d'xVh (RijW - ii? 2 ) , (3.15) 

where Rij is the Ricci tensor of the boundary metric and boundary indices i,j are raised and 
lowered with the boundary metric. This term is required to cancel logarithmic divergences as 
I, G — > 0. The finite part of this term is arbitrary, which is why we have included the constant 
(3. The integrand of this term is a combination of the Euler density and the square of the 
Weyl tensor. The former just contributes a constant term to the action but the latter may 
have important physical effects so we shall include it. For a pure gravity theory, adding a Weyl 
squared term to the action results in spin-2 ghosts in flat space but we shall see that this is not 
the case when the Yang-Mills theory is also included. The final counterterm Set is the finite R 2 
counterterm defined in equation 

When the four sphere boundary is unperturbed, the metric in the AdS region is 

ds 2 = V'idy 2 + sinh 2 y^dx' 'dx^), (3.16) 



4 In the prefactor of this equation, R refers to the radius of the sphere. In the integrand it refers to the Ricci 
scalar. 
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and the sphere is at y = y , where y is given by R = Isinh y . Note that y — ■> oo as I — > since 
i? is fixed. In order to use AdS/CFT for the perturbed sphere, we need to know how the metric 
perturbation extends into the bulk. This is done by solving the Einstein equations linearized 
about the AdS background. 

Our first task is therefore to solve the Einstein equations in the bulk to find the bulk metric 
perturbation that approaches on the boundary. We shall impose the boundary condition that 
the metric perturbation be regular throughout the AdS region. The most general perturbation 
of the bulk metric can be written 

ds 2 = P(dy 2 + sinh 2 y^ ij dx i dx i ) + Ady 2 + 2B i dydx i + H ij dx i dx j . (3.17) 

The first step is to decompose the bulk metric fluctuation into scalar, vector and tensor parts 
with respect to the four sphere: 

H i:j (y,x) = 9ij(y,x) + 2V^Xj)(y, x) + ViVj0(y,x) +'y ij ip(y,x). (3.18) 

The connection on the four sphere is denoted V. % is a transverse traceless symmetric tensor 
with respect to the four sphere: 

Vi0 ij ' = 6\ = 0, (3.19) 
where indices i,j are raised and lowered with 7^. \% is a transverse vector: 

V iX * = 0. (3.20) 

We can also decompose Bi into a transverse vector and a scalar: 

B t = B t + diB. (3.21) 

The quantities that we have introduced are gauge dependent. If we perform an infinitesimal 
change of coordinate then the five dimensional metric perturbation undergoes the gauge trans- 
formation 

Sg^ -> Sg^ + V M & + V„£ M - (3.22) 

We are using Greek letters to denote five dimensional indices. V is the connection with respect 
to the background AdS metric. The gauge parameters ^ can be decomposed with respect to 
the four sphere. £ y is a scalar and can be decomposed into a transverse vector and a scalar. 
Thus in total, we have four scalar degrees of freedom in our metric perturbation but there are 
two scalar gauge degrees of freedom so we can only expect two gauge invariant scalars. Similarly 
we have two vectors in our metric perturbation, but one vector gauge degree of freedom so there 
is only one gauge invariant vector quantity. The tensor part of the metric perturbation is gauge 
invariant. It is easy to check that the following scalar quantities are gauge invariant: 

1 y ( cosh y sinh y ) ' ^ ^ 

v& 2 = B - \d y <^ - / . - + cothy0. (3.24) 

2 2 cosh y sinh y 

Note that the residual gauge invariance discussed in section 3A. is also present here - we shall 
have more to say about this later on. 
The gauge invariant vector quantity is 

Xi = Bi- d y Xi + 2 cothyxi- (3.25) 

The gauge invariant tensor is %. 
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3.3 Solving the Einstein equations: scalars and vectors 

The Einstein equation in the bulk is 

1 6 

Rfiu - 2 R 9vv = p9nv- (3.26) 

We want to solve this such that our metric matches onto the perturbed metric on the four sphere 
boundary. The solution for the unperturbed sphere is simply AdS. Denote this background 
metric by g^ v . Linearizing around this background yields the equation 

V M V^ + V v V p 8g plt - V 2 8 9lMU - V^VJg? = ^8g, iu - ^g^8g p p , (3.27) 

This equation is gauge invariant and can therefore be expressed in terms of the gauge invariant 
variables. The yy component gives 

V 2 ^i - 2dyV 2 ^ 2 - 4 cosh y sinh yd y ^ x - 8 sinh 2 y\ff 1 = 0. (3.28) 

The vector part of the iy components gives 

V 2 X, = -3Xi. (3.29) 

The scalar part of the iy components gives 

di (cosh y sinh y^ x - 2^ 2 ) = 0. (3.30) 

The tensor part of the ij components gives 

dpij - 4 coth 2 yOij + cosech 2 ?/V 2 % = 0. (3.31) 

The vector part of the ij components gives 

(d y + 2cothy)V {i X j) = 0. (3.32) 

The scalar part of the ij components gives 

V i V i (-*i + 2d y * 2 + 4cothy# 2 ) 
+ % (cosh y sinh yd y ^ x + (8 cosh 2 ?/ - 2)*i + 2 cothyV 2 ^ 2 ) = 0. (3.33) 

Solving equation |3.32| yields 

V(iXj)(y,x) = Sin ^ 2 yo V(,Xj)(y ,g), (3.34) 
smh y 

which is singular at y — 0. We must therefore take the solution 

V {l X j) (y,x) = 0. (3.35) 
Thus the gauge invariant vector perturbation vanishes: we are free to choose X{ = 0. 
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Rearranging the equations for the scalars, one obtains 



V 2 ^i = -4*1 (3.36) 

and 

(coshysinhydy + (4 cosh 2 y - 2)) (ViVj + *i = 0. (3.37) 
This has the solution 

(V,V, • %) ^(y,x) = SinI ff CQS ^° (V,V i + %) Mvo,x). (3.38) 
v 7 smh ycosh y K ' 

Once again, this is singular at y = unless we take 

(V < V J -+7ii)*i(!/ > a;) = 0. (3.39) 

There is a regular solution to this equation, however it is simply an artifact of the ambiguity in 
our metric decomposition discussed in section [D] (see equation |3~5| ) so *i can be consistently 
set to zero. Equation |3.30| then implies that *2 is an arbitrary function of y. This is again 



related to an ambiguity in the metric decomposition: we are free to add an arbitrary function 
of y to <p without changing the metric perturbation. Hence we can choose *2 — 0. 

To summarize: we have solved the bulk Einstein equation for the gauge invariant vector and 
scalars, obtaining the result 

m 1 = ^ 2 = Xi = 0. (3.40) 

So far we have been working in a general gauge. We shall now specialize to Gaussian normal 
coordinates, in which we define ly to be the geodesic distance from some origin in our ball 
of perturbed AdS, and then introduce coordinates x % on surfaces of constant y (which have 
spherical topology). In these coordinates we have 

A = B = Bi = 0. (3.41) 

The presence of a metric perturbation implies that the boundary of the ball is not at constant 
geodesic distance from the origin. Instead it will be at a position 

y = ya + £{x). (3.42) 



We can now use our solution |3.40| to write down the bulk metric perturbation in Gaussian 
normal coordinates: 

ip(y, x) = f(x) sinhycosh?/, (3.43) 

(j)(y, x) = f(x) sinh y cosh y + g{x) sinh 2 y, (3.44) 

Xi(y, x) = Xi(x) sinh 2 y, (3.45) 

where /, g are arbitrary functions of x and x% ls an arbitrary transverse vector function of x. 
We now appear to have three independent scalar functions of x to deal with (namely /, g and 
£). These should be specified by demanding that the bulk metric perturbation match onto the 
boundary metric perturbation. However the boundary metric perturbation is specified by only 
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two scalars. We therefore need another boundary condition: regularity at the origin. Solutions 
proportional to sinh y cosh y are unacceptable since they lead to 



g^Sgau oc cothy, 



(3.46) 



which is singular at y = 0. We must therefore set f(x) = 0. To first order, the induced metric 
perturbation on the boundary is 



hij(x) = H i:j (y ,x) + 2/ 2 sinh y cosh y %£. 



(3.47) 



Recall that is given by equation |3.18| . The left hand side is decomposed into scalar, vec- 
tor and tensor pieces in equation [3.2p| . We can substitute the solution for the bulk metric 
perturbation into the right hand side and read off 



(3.48) 

(3.49) 
(3.50) 

These equations determine £(x), g{x) and Xi( x ) in terms of the boundary metric perturbation. 
In section [Oj , we showed that (f>(x) and Xi( x ) could be gauged away so we shall now set 



ip(x) = 2l 2 sinh y cosh y £ (x) 

(f)(x) = g(x) sinh 2 ?/o, 
Xi(y, x ) = Xi( x ) sinh 2 ?/. 



9{ x ) 



0. 



Xi{ x ) 



0. 



This implies that 



<f){y,x) = ip(y,x) 



0. 



o. 



(3.51) 



(3.52) 



In other words, all scalar and vector perturbations vanish in the bulk: the bulk perturbation is 
transverse and traceless. The only degrees of freedom that remain are therefore the bulk tensor 
perturbation and the scalar perturbation £(x) describing the displacement of the boundary. 



3.4 Tensor perturbations 

The tensor perturbations are less trivial: we have to solve equation |3.31|. This was done in [E4 



by expanding in tensor spherical harmonics These obey 

0, 

and they are regular tensor eigenf unctions of the Laplacian: 



f j H%\x) 



w i m'\x) 



V 2 H^ = (2-p(p + 3))H%\ 



(3.53) 



(3.54) 



where p = 2,3,.... We have suppressed extra labels k,l,rn,... on these harmonics. The 
harmonics are orthonormal with respect to the obvious inner product. Further properties are 
given in 



28 



5 We apologize for our slightly confusing notation: ip( x )> 4>{ x ) an d Xi( x ) i n equation |3.2| have, so far, nothing 
to do with the bulk quantities ip(y,x), 4>(y,x) and Xi{lJi x )- 
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The boundary condition at y = y isf] 9ij(y ,x) = where %(x) is the tensor part of 

the metric perturbation on the boundary. Imposing this condition together with regularity at 
the origin gives a unique bulk solution |24j 

%M = E^^W / dVyVV^V), (3-55) 
where / p is given in terms of a hypergeometric function: 

fp(y) = or Z 2 i 71 i(p/2,(p+l)/2,p + 5/2,tanh 2 y). (3.56) 

cosrr y 

3.5 The gravitational action 

We have now solved the Einstein equations in the bulk and found a solution that matches onto 
the metric perturbation of the boundary. The next step is to compute the action of this solution. 
The bulk contribution from the Einstein-Hilbert action with cosmological constant is 

I 3 C * I [Vo+£ 

Sbuik = J d xyj dy sinh y (3.57) 

- J d 5 xy/g~ - (r^ - ^Rg^u - p9ia^j Sgf^ - Sg^A^Sgpa 

The term that is first order in 5g^ u will vanish because the background obeys the Einstein 
equation. The second order term involves the Lichnerowicz operator (generalized to include the 
effect of a cosmological constant) Al, which is a second order differential operator with the 
symmetry property 

Af pCT = Af (3.58) 
This term vanishes because the perturbation is on shell, i.e., 

Af* 7 ^ = 0. (3.59) 

We are left simply with the background contribution 

Sbuik = TTTy f d A xJ^f [ dy sinh 4 y (3.60) 

Z7TG J JO 

l 3 £l d fvo , . . I 



2nG 



pyo / p i~~ 

dy sinh 4 y + J d 4 x^J 7 (4 sinh 4 y £ + 8 sinh 3 y cosh y £ 2 ) 



where denotes the volume of a unit four sphere. Of course, in order to rearrange the 
Einstein-Hilbert action into the form |3.57 we have to integrate by parts several times, giving 



rise to surface terms. These will depend on derivatives of the bulk metric perturbation evaluated 
at the boundary. Since there are only tensor degrees of freedom excited in the bulk, only tensors 
will occur in these surface terms - there will be no dependence on £. The surface terms are 



Ssurf = J^J ( y^d y B l3 - ^2^0% ) . (3.61) 



6 The boundary is actually at y = yo + which gives higher order corrections. These would appear at 
third order in the action as couplings between tensors and scalars. 
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The second contribution to the gravitational action is the Gibbons-Hawking term. In 
evaluating this, it is important to remember that the unit normal to the boundary changes 
when we perturb the bulk metric. The boundary is a hypersurface defined by the condition 
f(y,x) = y — £(x) = yo- The unit normal is therefore given to second order by 

n = l{ 1 - ^ff ) dy - Td£dx\ (3.62) 
\ 2 sinh y J 

Note that this holds for a range of y and therefore defines a unit covector field that is normal 
to the family of hypersurfaces / = constant. In other words, it defines an extension of the unit 
normal on the boundary into a neighbourhood of the boundary. Written as a vector, the normal 
takes the form 

n = \L_ d&H \ d_ _ ( H _ e^(y,x)d^ \ _d_ 
I \ 2 sinh 2 ?// dy \lsmh 2 y I 3 sinh 4 y / dx v 

where #y (y, x) is the bulk tensor perturbation. The trace of the extrinsic curvature is 

K = V/. (3.64) 

In evaluating this one must take account of both the perturbation in the unit normal and the 
perturbation in the connection. The result is 



K = Uothy- T ^n-^-d i m 
I I sinh y I sinh y 



This has to be evaluated at y = y + ^. To evaluate ^fy on the boundary, we need to know the 
induced boundary metric perturbation to second order: 

h i:j (x) = 6ij(y , x) + 2P sinhyo coshy 7ij^ + P (2 sinh 2 y + l) + Pd&djt + £d y 9 i:j . (3.66) 
These results can now be substituted into the Gibbons-Hawking term, yielding 
P 



Sgh = ~~^q J <^ x \jl 4 cosh y sinh 3 y + sinh 2 ?/ (l6sinh 2 ?/o + 12) £ (3.67) 
+ cosh y sinh y (32 sinh 2 y + 12^ £ 2 - 3 cosh y sinh y £V 2 £ - ^0 %j d y 8ij . 

We have integrated some terms by parts. So far, we have expressed the scalar part of the action 
in terms of £. However, we really want to express everything in terms of the induced metric on 
the boundary, which has scalar part ip(x). This can be done by taking the trace of equation 
3.66| and solving for £ in terms of ip to second order, giving 



(2sinh^ + ljr dd&il) 



2l 2 sinh y cosh y 8Z 4 sinh y cosh y 32Z 4 sinh y cosh y 
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The total contribution from the Einstein-Hilbert and Gibbons-Hawking terms is given by the 
sum of the following 



o(o) 

u grav 



3m 4 



2nG Jo 



dy sinh 2 y cosh 2 y, 



Q(l) 
grav 



3/ 3 



AttG 



J d A x^p{= cosh yo sinh y ^, 



(3.69) 
(3.70) 



C(2) 
u grav 



8nG 



d x\/7 



3 2sinh z y + 1) ip 



3^V 2 ^ 



2Z 4 sinh y cosh y 8/ 4 sinh y cosh y 



8/ 4 * 1 2/ 4 



'./ 



(3.71) 



We can now expand the action in powers of l/R (using sinhy — R/l)- This gives terms that 
diverge as /~ 4 and l~ 2 as J goes to zero. For the scalar perturbation, these divergences are 
cancelled by the counter terms Si and S2. For the tensor perturbation (dealt with in [p^fl), the 
third counter term S3 is needed to cancel a logarithmic divergencef]. 

The final term that we have to include in the effective action is the finite counter term Set- 
Evaluating this to second order gives 



s 



ct 



3aN 2 Q 4 3aN 2 
+ 



4tt 2 64tt 2 j R 4 
The final result for the Yang-Mills effective action is 

w = w {0) + w {1) + w {2) + ... 

where 



J d 4 x^fi Uv^ + 4^v 2 ^ + -e ij v 2 e lj - g#%) • (3.72) 



w {0) 



+ -^r^ + ~^r^ (41og2 - 1) 



8tt 2 



Att 2 
3N 2 
16n 2 R 2 



32tt 2 



(3.73) 

(3.74) 
(3.75) 



where 



3N 2 



+ 



Q4n 2 R 4 
N 2 

256n 2 R 4 



J d A x^i [V> (V 2 + 2) i) - aip (V 4 + 4V 2 ) ip 



x (*(p) + 2/3p(p + l)(p + 2)(p + 3) - 4ap(p + 3)) , 



*(P) = p(j9+l)(p + 2)(j9 + 3)[^(j9/2 + 5/2)+^(p/2 + 2)-^(2)-^(l)] 
+p 4 + 2p 3 - 5p 2 - 10p - 6. 



(3.76) 



(3.77) 



7 This counter term is formed from the Euler number and the square of the Weyl tensor, neither of which is 
affected by scalar perturbations. S3 therefore does not contribute to the action for scalar perturbations. 
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The scalar perturbations have an action that can be expressed simply in position space. How- 
ever, the tensor perturbations are given in momentum space where they have an action with 
complicated non-polynomial dependence on p. This corresponds to a non-local action in position 
space. At large p it behaves like p 4 logp, as expected from the flat space result for (Tij(x)Ti>j>(x')) 
0. 



3.6 Metric correlation functions 

Our theory is just four dimensional Einstein gravity coupled to the Yang-Mills theory, with 
action 



S 



±^Jd 4 x^R+W, (3.78) 
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where we have not included a Gibbons-Hawking term because the instanton has no boundary. 
Note that we are still working in Euclidean signature. W denotes the Yang-Mills effective action, 
including the effect of the finite counterterms. G is the four dimensional Newton constant. In 
order to compute the two point correlation functions of metric perturbations we need to calculate 
the terms in S that are quadratic in the metric perturbations described by and ip. 
To second order, the Einstein-Hilbert action of the perturbed four sphere is 

(3.79) 

Adding the Yang-Mills effective actions gives the total action. This has a non-vanishing piece 



linear in ip. Varying ip fixes R to take the value given by equation 2.19, which implies that the 



linear term vanishes. Equation 2T9 can be used to write G in terms of R, which brings the 



quadratic part of the scalar action to the formQ 

S scalar = J ^1 4 faV 2 ~ l) (V* + 4) ^ } (3.80) 

and the quadratic part of the tensor action becomes 

N 2 _ , ' " - ^ - 

V 

where 

F(p, a,p)=p 2 + 3p + 6 + ^(p) + 2{3p(p + l)(p + 2)(p + 3) - Aap(p + 3). (3.82) 
From these expressions we can read off the correlation functions of metric perturbations: 

1 1 



Stensor = ^E(/ d*x' ^ 9* (x')!!^ (x'tf F(p, a, (3) , (3.81) 



(tp(x)?p(x')) 



3N 2 (-a)(A + m 2 



m? 



(3.83) 



8 If a = then this is almost exactly the same as the scalar action one would obtain for perturbations about 
a de Sitter solution supported by a cosmological constant. The only difference is that the overall sign is reversed. 
This implies that, with the exception of the homogeneous mode, the conformal factor problem of Euclidean 
quantum gravity is solved by coupling to the Yang-Mills theory when a — 0. 
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where 



The tensor correlator is 



m 2 = — . (3.84) 
2a K ' 



(%(x)M*')> ^E^j'l 1 .*^)" 1 . ( 3 - 85 ) 

where the bitensor W]jht(x,x') is defined as 

W$ f (x,x>)= Y, H%\x)H$(x'), (3.86) 

k,l,m,... 

with the sum running over all the suppressed labels k,l,m, . . . of the tensor harmonics on the 
four sphere. 

4 Analytic structure of propagators 
4.1 Flat space limit 

Before analyzing our correlation functions we shall consider the analagous functions in flat 
space. This will allow us to constrain the allowed values of the parameters a and (3, which will 
be important when we return to the de Sitter case. 

Recall that in equations |3.75| , |3.76| and |3.79| , the radius R is arbitrary. To avoid confusion, 



we shall now denote this arbitrary radius by R to distinguish it from the on-shell value R, given 
by equation |2.19| . We can recover flat space results by taking R — > oo. Before taking this limit, 
we first replace the dimensionless momentum p with the dimensionful momentum k = p/R. 

There is no conformal anomaly in flat space and the scalar ip corresponds to a conformal 
transformation. Therefore, the only matter contribution to the scalar propagator comes from 
the term in the Yang-Mills action that breaks the conformal invariance, namely the finite counter 
term S ct . The other contribution to the scalar correlator comes from the Einstein-Hilbert action. 
One obtains 

M*M*0> oc _ Q2 \ M2 - (4.1) 
with a positive constant of proportionality. M 2 is given by 



where R is given by equation [2.19| , although we emphasize that we are now working in flat 
space. The second term in the propagator describes a massless scalar ghost. This can be dealt 
with by gauge fixing the action. The first term is more worrying. If a > then it describes a 
tachyon. We regard this as undesirable: we do not want flat space to be an unstable solution 
of our theory. We shall therefore always take a < 0, which gives a massive scalar in flat space. 

For the tensor propagator, the limit R — > oo makes the coefficient of the third counterterm 
S3 diverge. To cancel this divergence, introduce a length scale p defined by 

P = \og(p/R). (4.3) 
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The R dependence in the coefficient of the third counterterm then drops out, leaving a finite 
coefficient depending on the renormalization scale p. The R — > oo limit of the propagator is 
similarly well-defined. The result is proportional to 

(4.4) 



k 2 {1 + R 2 k 2 [l + log(fcV/4)]}' 



Our propagator is of exactly the same form as given by Tomboulis in his analysis of the 
effects of large N matter on the flat space graviton propagator. The propagator is defined for 
k 2 > 0. It can be analytically continued into the complex k 2 plane by taking a branch cut 
for the logarithm along the negative real axis. There are generally two poles present, with 
positions dependent on p. If p < 2R/e then these poles are on the positive real axis. One has 
positive residue and the other negative residue, so they correspond to a tachyon and a ghost. 
As p — > 2R/e, the two poles move together and merge to form a double pole. For p > 2R/e, 
this double pole splits into a pair of complex conjugate poles which move off into the complex 
k 2 plane. The modulus r and phase 9 of k 2 at these poles are related by 

r = w (4 ' 5) 

9 is given by solving 

Co 2 sin 9 \ 

1 + log^ + log— J, (4.6) 

which is straightforward to analyze graphically. The solution obeys 9 — > ±7r and r — > as 
p — > oo. 

The presence of tachyons for small p was not mentioned by Tomboulis since he implicitly 
assumed p ^> R. Since we want flat space to be a stable solution of our theory, we shall take 
p > 2R/e when we consider the propagator in de Sitter space. This corresponds to taking 
> log 2- 1. 

It is interesting to note that changing p changes the coefficient of the third counter term 
S3 by a finite amount. This corresponds to introducing a finite counter term involving the 
Euler number and the square of the Weyl tensor. The former is left unchanged by metric 
perturbations. However, the latter is known to give rise to spin-2 ghosts in a pure gravity 
theory. Such ghosts do not appear in our model: coupling to the CFT removes them. 

4.2 Scalar propagator on the sphere 



Equation |3 . 8 3| is the propagator of scalar metric perturbations on a spherical instanton supported 



by the conformal anomaly of the CFT. The first term in the propagator describes a particle 
with physical mass-squared m 2 /R 2 = (2aR 2 )~ 1 . Since we are assuming a < 0, we have m 2 < 
so this particle is a tachyon. This is good because we do not want the spherical solution to be 
stable since that would lead to a Lorentzian de Sitter solution in which inflation never ends. 
Making a more negative makes the tachyon mass squared less negative, and therefore makes 
the instability weaker. This suggests that if a is sufficiently negative then inflation will last for 
a long time. We shall make this more precise later. 

The second term in the propagator describes a ghost. This is the normal scalar mode of 



gravity that is canceled by the scalar parts of the Fadeev- Popov ghosts [ 26| . These ghosts supply 
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a determinant that cancels the (V 2 + 4) factor in the scalar action. The propagator can then 
be read off from the action: 



(1>(*)4>W) = (-V 2 + m 2 )" 1 . (4.7) 

This propagator can be written in momentum space as 

;/ 3|a|iV 2 ^ p(p + 3) + m 2 ' V ; 

where the biscalar is a function of the geodesic distance fi between x and x', given by 

W (p) (/i(x,x')) = E H^{x)H (p \x'), (4.9) 

k,l,m 

where if ( p ) denote spherical harmonics on the four sphere and the sum runs over the suppressed 
eigenvalues k, I, m. 

Notice that there are many negative modes if a is negative and close to zero. However, 
if a < —1/8 then only the homogenous (p = 0) negative mode remains. To compute the 
primordial density fluctuations in the microwave background radiation we are interested in the 
two-point function with the homogenous mode projected out [3~7|. Notice also that the Fadeev- 
Popov ghosts fix the residual gauge ambiguity associated with the p — 1 modes. These modes 
no longer have zero action and therefore cannot be regarded as gauge. 



4.3 Tensor propagator on the sphere 

The tensor propagator (equation |3.85|) has an interesting analytic structure. The momentum 
space propagator is proportional to F(p, a, /3) _1 , where F is given by equation \5.82\ 

For a physical interpretation, we need to study the behaviour of F in the complex X p plane, 
where X p = p(p + 3) — 2 is the eigenvalue of —V 2 . We must therefore first write the propagator 
as a function of X p . Since 

P=-l±^ + X P , (4.10) 

we must choose a branch for the square root. The Euclidean propagator is defined as a sum 
over p = 2, 3, . . ., for which X p is positive. We must therefore take the positive sign for the 
square root. The analytic continuation into the complex X p plane is given by taking a branch 
cut along the negative axis for X p < —17/4. p has positive imaginary part just above the cut 
and negative imaginary part just below the cut. Note that Re(p) > —3/2. The branch cut 
corresponds to a continuum of multi-particle states. The imaginary part of the propagator is 
discontinuous across the cut. In general, the absence of negative norm states implies that the 
imaginary part of the propagator just below the cut minus the imaginary part just above the 
cut should be positive, which is indeed the case for our tensor propagator. 

It is also possible for the tensor propagator to have discrete poles in the A p plane. Poles on the 
real axis are of particular importance. If such a pole occurs at positive X p then it corresponds 
to a tachyon. In fact, since the graviton in de Sitter space has an equation of motion with 
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Figure 2: Inverse propagator F(p, 0,0) for —3/2 < p < 1/2 and —0.1 < p < 0.1. The graph 
grows monotonically for p > 0. There are zeroes at p ~ —1.48 (massive particle), p ~ —0.054 
(ghost) and p = (massless graviton). 



A p = —2, it seems appropriate to regard particles with X p > —2 as tachyons. If a pole on the 
real axis has negative residue then it corresponds to a ghost. 

Our propagator always has a pole at \ p — — 2 (p = 0), corresponding to the massless 
graviton in de Sitter space. Support for this interpretation comes from observing that transverse 
traceless tensor harmonics have 5 degrees of freedom. However, the mode with p = mixes 
with transverse vector harmonics, which have 3 degrees of freedom. Thus the p = mode has 
3 gauge degrees of freedom, leaving 2 physical degrees of freedom, as appropriate for a massless 
spin-2 particle. 

We shall start by considering the case a — (3 — 0, for which there are two other poles in our 
propagator, one at p « —1.48 and the other at p « —0.054. The former has \ ~ —17/4 (but is 
not quite on the cut) and has positive residue, the latter has \ ~ —2.16 and negative residue. 
The behaviour of F(p, 0, 0) is plotted in figure It is easy to show that signs of the residues 
of F _1 with respect to \ p are given by the slope of F as it passes through 0. The positions of 
the poles are shown in figure ^. Changing the value of (3 (still with a = 0) changes the position 
and nature of these poles. As (3 is made more positive, the pole with p « —1.48 gets absorbed 
into the branch cut and the ghost moves towards p = — 1 (i.e. X p = —4). As (3 is made more 
negative, the pole with p « —1.48 moves towards p — — 1 while the other pole moves to positive 
p (i.e., A p > —2), with its residue changing sign as it crosses p — 0. This pole corresponds to a 
tachyon. Recall that tachyons were also present in flat space for sufficiently negative j3. In order 
for tachyons to be absent in flat space, we had to choose (3 > log 2 — 1. We have roughly the 
same restriction on (3 in order to avoid spin-2 tachyons in de Sitter space. We shall therefore 
exclude the case (3 < log 2 — 1 as unphysical. 

Now consider the effect of turning on a < 0. This has no effect on the pole at A p = —2, so 
the massless graviton remains. If (3 — 0, then the two other poles move together as a decreases 
and eventually coalesce into a double pole. This splits into a pair of complex conjugate poles 
that move off into the complex A p plane. For (3 > then there is generally only one pole present 
(in addition to the graviton pole) when a = 0. As a is decreased, an additional pole (with 
positive residue) emerges from the branch point and moves towards the ghost pole, eventually 
coalescing with it. This then splits into a pair of complex conjugate poles. If (3 < 0, then the 
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Figure 3: Analytic structure of the tensor propagator in the complex X p plane when a = 0. 
The dotted lines denote X p = —2. Poles on the real axis to the right of this line correspond to 
tachyons. There is a branch cut at X p = —17/4 and the thick line represents the branch cut. 
There is always a massless graviton pole at X p = —2. The diagram on the left is for (3 > 0, 
when there is a single ghost pole. As (3 decreases, this pole moves to the right and another pole 
emerges from the branch cut. This new pole corresponds to a massive particle and appears in 
the second diagram, which is for (3 = 0. The final diagram is for (3 < 0, when the ghost pole 
crosses through X p = —2 and becomes a tachyon. 



two poles again move together, coalesce and then become a pair of complex conjugate poles. 
In all cases, the effect of making a more negative is similar to the effect of increasing p in the 
flat space propagator, i.e., pathologies such as ghosts and tachyons move off into the complex 
plane. When (3 is large, the poles becomes complex for a < —(3/8, so no fine tuning of the ratio 
at/ (3 is involved. 



4.4 Complex poles 

We have seen how ghost poles can be moved off the real axis, becoming a pair of complex 
conjugate poles. The interpretation of such a pair of poles has been reviewed by Coleman |38[ . 
The presence of complex conjugate poles with (complex) masses given by m = a ± ib with 
b > implies causality violation at lengths or times of the order of 1/ V&. For Tomboulis' flat 
space propagator, we have b ~ -R -1 , so one expects causality to be violated at a length scale 
of the order of R, which is roughly iV times the Planck length. Unless N is enormous, this is 
far less than any scale probed by particle physics experiments so such causality violations are 
unobservableQ, as noted by Tomboulis. 

For our de Sitter propagator, the complex poles again have b oc R" 1 . If \a\ is large then b oc 
\/— <xR _1 , so causality violation occurs on a time scale Rj \f—a. If \a\ is not large then causality 
violation occurs on a time scale R. This is much smaller than scales probed in experiments, 
but may have observational consequences in the CMB since R is the Hubble time, and therefore 
the time scale for microphysics during inflation. However, we shall see in the next section that 
observations suggest that |a| is of order 10 9 , so causality violation occurs on a time scale much 
shorter than the Hubble time and is therefore completely unobservable. 



9 In fact, these effects might be smaller than the effects of the gravitational field of subatomic particles, which 
would also lead to modifications of causality through tilting of light cones. 
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5 Lorentzian two-point correlators 



In this section we will show how the scalar and tensor propagators on the four sphere instanton 
uniquely determine the primordial CMB perturbation spectrum in Lorentzian closed de Sit- 
ter space. The two-point correlators in the Lorentzian region are obtained directly from the 
Euclidean propagators by analytic continuation. We refer the reader for the details of this cal- 
culation to our previous paper |24|] , where we described the analytic continuation of the graviton 
correlator in a Randall- Sundrum version of the Starobinsky model. The techniques to perform 



these calculations were developed in [39, 40 



5.1 Scalar propagator 



We have the Euclidean correlator |4.8| as an infinite sum over real p, where p labels the level of 
the four sphere scalar harmonics. Although this is a convenient labelling to study their analytic 
structure, the eigenspace of the Laplacian on de Sitter space suggests that the Lorentzian 
propagator is most naturally expressed in terms of an integral over real positive p' = i(p + 3/2), 
corresponding to scalar harmonics of the Lorentzian Laplacian with eigenvalue X p i = (p' 2 + 9/4). 
We must therefore first analytically continue our result for the propagators into the complex 
p-plane before continuing to Lorentzian signature. In terms of the label p', the Euclidean scalar 



correlator 4.S becomes 



with 



m m " 3|a|iV 2 p; ^ /2 p' 2 + 9/4-m 2 ' 1 ; 

W V) {z) = 5 «PV + V4) 2Fi(3/2 + 3/2 _ ^ g) (5 2) 

and z = cos 2 (/i/2). This biscalar is analytic in the upper half p'-plane. The coefficient of the 
biscalar is also analytic in the upper half plane apart from a simple pole at p' = A t , where 



A t = iJ--m 2 . (5.3) 



This pole corresponds to the tachyon. Notice that the sum in equation 5T starts at p' = hi/ 2 
because we have projected out the negative homogenous mode, which should be regarded as 
part of the background f37| . 

Knowing the analytic structure of the correlator, we are able to write the sum [5J] as an 
integral along a contour C\ encircling the points p' = 5z/2, 7i/2, ..to/2, where n tends to infinity. 
This yields 

Wfl) * ( " }> = k V + 9/4 - ^ • (5 ' 4) 

The contour C\ can be distorted to run along the real p'-axis. Apart from the tachyon pole, 
we encounter two extra poles at p' = 3i/2 and p' = i/2 in the tanhp'7r factor. The p' = 3i/2 
pole corresponds to the negative homogenous mode that we have projected out in the Euclidean 
correlator. On the other hand, W^^(fjL) = so the pole at p 1 = i/2 does not contribute to 
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the propagator. The contribution from the closing of the contour in the upper half p'-plane 
vanishes. Hence our final result for the Euclidean correlator reads 

16i7r 2 R 4 



3\a\N 2 



+°° , ,{tanhi/ir)WW(z) 
dp 



p' 2 + 9/4 - m 2 



a7 



(tanhA f vr)H/ (At) (^) + 



lOz 



m 2 7r 2 



(5.5) 



Finally one can rewrite |5^ as an integral from to oo, over the eigenspace of the Lorentzian 
Laplacian, and the two discrete contributions from the tachyon pole and the homogenous mode. 
The tachyon contribution grows exponentially for timelike intervals. However, the relevant 
propagator for computing the CMB anisotropies is the Feynman propagator, which should be 
bounded both to the past and future. Therefore, the propagator that we have obtained by 
analytic continuation from the four sphere does not obey the appropriate boundary conditions. 
In order to obtain the two-point function that describes the correlations in the primordial density 
fluctuation spectrum, we change the contour of integration so as to exclude the contribution 
from the tachyon pole. We then obtain the Lorentzian Feynman scalar propagator, 



(ijj(x)ip(x')) 



32n 2 R 4 
'3|a|iV 2 



o 



+ °° dp' ( tann P' 7r ) W/L(p ' ) (^( a; ) x ')) 



p' 2 + 9/4 



+ 



10 



vnr 



m 2 ir 2 



(5.6) 



The Lorentzian biscalar W L ( p/S) differs from W^' only by a factor of —i and (tanh.p , 7r)W L ^ p '\z) 
equals the sum of the degenerate scalar harmonics on closed de Sitter space with eigenvalue 
X p i = {p' 2 + 9/4) of the Laplacian. For spacelike separations, we have z = cos 2 (/z/2), where 
fj,(x, x') is the geodesic distance between x and x' . The correlator for timelike intervals is 
obtained by setting p = ir/2 — it, where p is the polar angle on the four sphere. For a purely 
timelike separation, this gives z = cosh 2 ((t — if)/ 2). 



5.2 Tensor propagator 



The principles of the continuation of the tensor propagator [3.85| are the same, but the calculation 



is more complicated. We refer the interested reader to our previous paper [24| for the technical 



details. The differences between [^3] and the present paper are that we now have included the 
effect of the finite R 2 counterterm, we have kept (3 in the coefficient of the third counterterm 
arbitrary and we now treat the discrete poles in the propagator more carefully. 

In [24]] it was shown that the bitensor W^2,(/i) can be unambiguously extended as an 



analytic function into the upper half p'-plane. In addition, from subsection ^]3| we know that its 
coefficient F(—ip' — 3/2, a, /3) _1 is analytic, apart from a simple pole at p' = Si/2, corresponding 
to the massless graviton in de Sitter space, and a pair of poles with complex masses Ai and 
A 2 = — Ai (we are assuming that a < — /3/8 so that there are complex poles instead of a ghost). 
These poles always occur in the upper half p'-plane. 



Writing the sum in equation |3.85| as a contour integral yields 

6Mir 2 R 4 r . , t{v >) 



N 2 



J c dp'ttmhp'irW^jr&Gip', a, /?)" 



(5.7) 



where 



G(p',a,(3) = F(-ip' -3/2, a, 0) 

= p' A - Aip' 3 - p' 2 /2 - hip' - 3/16 + (p 12 + 9/4) [4a + (p' 2 + 1/4) x 
(ip(-ip'/2 + 5/4) + ij){-ip'/2 + 7/4) - ^(1) - ^(2) + 2/5)]. 
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As we deform the contour towards the real axis we encounter, apart from the poles mentioned 
above, two extra poles in the tanhpV factor. However, as explained in detail in ||24j| , they do not 



contribute to the tensor fluctuation spectrum. The contribution from the closing of the contour 
in the upper half p'-plane vanishes. Using G(—p', a, (3) = G(p', a, 0), one can again rewrite the 
remaining integral over the real axis as an integral from to oo. The continuation of z(x,x') 
for timelike intervals is the same as for the scalar two-point function. We then obtain for the 
Lorentzian tensor propagator, 

(MaOM**)) = -^-{l ^ , (tanhpV)^j,(^)K(G(p',a,/?)- 1 ) 



irRiji'fiz) - 2vrK [(tanh A^W^iz) R (Al) ] } • (5.8) 



In the integral, (t&nhp'ir)W L Mj,(z(x,x')) can be identified with the sum of the degenerate 
rank-two tensor harmonics on closed de Sitter space with eigenvalue X p i = (p' 2 + 17/4) of the 
Laplacian. The integrand vanishes as p' — > 0, so the correlator is well-behaved in the infrared. 

The first term in equation |5.8j represents the continuous tensor fluctuation spectrum. The 
second term describes the massless graviton with Rj 3 yj-/(z) defined as the residue at p' = 3i/2 
of 

tanhpV 

" G(p>,a,f3y 

The third term in |5.8| is the combined contribution from the complex poles, with R(Aj) denoting 
the residue of G(p', a, (3)~ l at p' = A%. For large \a\ this mode grows exponentially, imply- 
ing that the analytically continued propagator does not obey the boundary conditions for the 
Feynman propagator. This can be remedied by changing the contour of integration to exclude 
the contribution from the complex poles, giving the correct propagator for two-point tensor 
correlations in the microwave background: 



WW A*) ru7i^ - ( 5 - 9 ) 



t-tji'j'yZ) 



1 9S7r 2 /? 4 r r+oo 

(%(x)^v(x')) = N2 J o dp'(tanhp'7r)W^j,(z) M(G(p', a, (3)' 1 ) - ttR, 

(5.10) 

If \a\ is large then the tensor propagator is proportional to (|a|iV 2 ) 1 . At large p' the tensor 
propagator behaves like (p' A logp') -1 , just as the Euclidean correlator |3.85| . This is in contrast 
to the usual p'~ 2 behavior of the graviton propagator for de Sitter space with a cosmological 
constant. 

6 Observational constraints 
6.1 Duration of inflation 

The Starobinsky instability in four dimensions has been analyzed carefully by Vilenkin |T7j| . He 
showed that the scale factor grows exponentially until 

6-fz n / \ , 

* = *.~^(7-l), (6-1) 
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where, for our model, the parameters H and M are given by 



H = R-\ M = {^/^2aRj . (6.2) 
The parameter 7 is related to the initial perturbation from the exact de Sitter solution 

7 =^log(2/5 ), (6.3) 



where 

, H -H 



(6.4) 



Ho ' 

is the perturbation of the Hubble parameter H = a/a at time t = 0. If 5o < then the solution 
eventually becomes singular jnj, at least if one neglects spatial curvature (which should be a 



good approximation if there is a lot of inflation). We shall therefore restrict ourselves to 5q > 0. 

For t < £*, there is exponential growth with Hubble parameter H . The number of e-foldings 
of inflation during this phase is therefore 

^1 = ^(7-1). (6-5) 

For our values of H and M, this gives 

Ni = -12a(7- 1). (6.6) 

For t > i*, there is a phase of slow- roll inflation in which the Hubble parameter changes 
from Hq to M. The number of e-foldings of inflation during this phase is JL7] 



iV 2 = -12a log cosh 1 » -2.26a. (6.7) 

The slow-roll phase lasts until t ~ Q^Hq/M 2 . Once this phase ends, the universe enters a matter 
dominated era in which the scale factor behaves as ITS, IT) 



a(t) oc t 2/3 1 + — — sin Mt + 0{r 2 ) . (6.8) 

The oscillations in the scale factor can drive particle production and reheating. 

Vilenkin used the Wheeler-DeWitt equation to obtain an estimate for 8 . Using his results, 
we obtain 

5 ~-i^, N 1 = -12a{logN-l). (6.9) 

Quantum cosmology therefore predicts 7>1. So far, the only restriction on iV is that iV must 
be large enough for our AdS / CFT calculation to be valid. This implies that log iV is not close 
to 1, so taking a < —5 makes Ni sufficiently large to solve the horizon and flatness problems. 

Our correlation functions for metric perturbations were calculated assuming a four sphere 
(or de Sitter) background. The present day horizon size left the horizon about fifty e-folds before 
the end of inflation. Hence the long-wavelength temperature fluctuations in the microwave sky 
carry the imprint of the first expansion phase provided N2 < 50, which is true if a > —20. 
Because our correlation functions for metric perturbations were calculated assuming a four 
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sphere (or de Sitter) background, the predicted spectrum can then be directly compared with 
observation. However, our results will be modified for modes that left the horizon during the 
slow-roll phase, when the background is not exactly de Sitter. Therefore, if a < —20 then it 
would be necessary to do a calculation based on an scalar/vector/tensor decomposition on the 
three sphere in order to enable us to evolve the spectrum through the instability and predict in 
detail the CMB fluctuation spectrum. 

6.2 Amplitude of perturbations 

In order to compare our results with observations, we should first render the propagators di- 
mensionless by dividing by R A . The correlators are then functions of p divided by iV 2 . Long 
wavelength perturbations are insensitive to what happens after inflation, so these can be di- 
rectly compared with observation. For the tensors, long wavelength perturbations correspond 
to modes on the four sphere^ with p = 2. The amplitude of the fluctuations can be obtained 
from the correlator: 

In order to agree with observations this should not exceed 10~ 5 , which requires 

iV 2 (250 + 240/5 - 40a) > 10 13 . (6.11) 

Since we are assuming N is large, the obvious way to satisfy this inequality is to take N = 
O(10 5 ). However, this implies that the number of fields present is lliV 2 = O(10 u ), which 
seems to contradict present day observations^. Instead, we could take N 2 f3 to be of order 
4 x 10 10 or iV 2 |a| to be of order 2 x 10 11 . The former corresponds to taking the coefficient of 
the Weyl squared term in the action to be of order 10 7 and the latter corresponds to taking the 
coefficient of the R 2 counterterm to be of order 10 8 . 

Note that if we take (3 to be large then we would also have to take a to be large in order 
to avoid ghosts in the tensor propagator. Therefore the most natural choice is probably to 
take just a to be large. Note that suppression of tensor perturbations through a Weyl squared 
counterterm (i.e. taking (3 large) was not mentioned in |15| [nj since this counterterm does not 
affect the coefficients a, c, d in the trace anomaly. 

Turning to the scalar perturbations, we see that these can also be suppressed by taking 
iV 2 |o;| to be large. Changing (3 does not affect the scalars. Our scalar correlator suggests that 
taking iV 2 |a| to be of order 2 x 10 11 should bring the scalar perturbations within observational 
bounds. 

We conclude that if iV 2 |o;| is of order 2 x 10 11 then we can bring metric perturbations within 
the observational bounds. iV just has to be large enough to justify the large N approximation 
for the matter fields. For example, we could take N = 10 and a = —2 x 10 9 . However, such a 
large value for a implies that all modes that we observe today must have left the horizon during 

10 We should really be studying the Lorentzian correlators here. However, the overall amplitude of the 
Lorentzian and Euclidean propagators is the same. 

11 However, it is possible that these fields may have masses large compared to the scale probed in colliders, 
i.e., m S> 1 TeV, but small compared with the scale at which inflation takes place, m <C 10~ 5 m p i. Such fields 
would be effectively massless during inflation but unobservable today. 
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the slow-roll phase of inflation. Our results for the two-point correlators will be modified in 
this case, since we assumed a four sphere background in our calculation. However, it is usually 
the case that the amplitude of perturbations is inversely proportional to the horizon radius at 
which they left the horizon. The horizon radius increases during slow-roll so it seems likely that 
if | a | is very large the amplitude of perturbations will be smaller than the amplitude obtained 
above. This argument is confirmed by the estimates of Vilenkin [I7|. We conclude that taking 



iV 2 |a| s;2x 10 11 will bring the perturbations within observational bounds, and a far smaller 
value may in fact be sufficient. 

A coefficient of order 10 8 in the action is large, but this is essentially the same fine-tuning 
problem that also appears in all scalar field models of inflation. In these scenarios, matching 
the amplitude of perturbations to COBE typically requires a fine-tuned parameter in the action 
ofO(10- 12 ). 

Note that taking \a\ to be very large implies that causality violations during inflation occur 
on a time scale much shorter than the Hubble time, so they would not have had a significant 
effect on microphysics. One might worry that taking \a\ to be large would imply significant 
deviations from Einstein gravity today, arising from the higher derivative R 2 term in the action. 
In flat space, the only effect of this term is to introduce a scalar field with mass given by equation 



4~2|. If we take N = 10 and led of order 10 9 then this scalar has mass M ~ 10 6 m„/, which is 



far too massive to be observed nowadays. 



7 Short distance physics 
7.1 Introduction 

The observational constraints that we have derived do not depend on the detailed structure 
of our propagators and could be obtained directly from the work of Starobinsky and Vilenkin. 
In this section we shall consider a new phenomenon revealed by our propagators, namely the 
suppression of short distance metric perturbations by matter fields. This suppression is evident 
in Tomboulis' flat space propagator |4.4j , which falls off as (fc 4 log fc 2 ) -1 for large momentum k. 
It is also present in our tensor propagator^] , equation [3.85| , which falls off as (p 4 logp) -1 at 
large p. This behaviour has not been discussed in previous studies of the Starobinsky model 
because these have neglected the non-local part of the matter effective action. 

Inflation acts as a "cosmic magnifying glass" by blowing up microscopic physics to macro- 
scopic scales. It is often assumed that this might lead to some characteristic signature in the 
CMB of new physics at short distances, e.g., extra dimensions. Our results appear to contradict 
this inflationary dogma, because they show that at small scales, matter fields will completely 
drown out the effects of any new gravitational physics. In this section we shall illustrate this 
phenomenon by comparing our results with the results for a model with an extra dimension, 
namely the Randall- Sundrum (RS) [^] version of the Starobinsky model. 



12 Once again, we shall concentrate on the Euclidean propagators in the section. The Lorentzian propagators 
exhibit similar short distance behaviour. 
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7.2 Randall-Sundrum model 



The RS model consists of a five dimensional spacetime with negative cosmological constant, and 
a thin positive tension domain wall whose tension is fine tuned to cancel the effect of the bulk 
cosmological constant. The ground state solution of this model is a Poincare symmetric domain 
wall separating two regions of AdS. In the RS version of the Starobinsky model, we simply add 
a U (N) Yang-Mills theory to the worldvolume of the domain wall. This model was extensively 
discussed in our previous paper |2~3 



For related work, see IO, §3L 03L P 
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where g av denotes the five dimensional bulk metric and the metric induced on the domain 
wall, I is the radius of the AdS solution, W is the generating functional of the Yang-Mills theory 
on the domain wall and [K]~t is the discontinuity in the trace of the extrinsic curvature at the 
domain walQ 

There are two simple solutions of the equations of motion for this model. Since the trace 
anomaly vanishes in flat space, a Poincare symmetric solution still exists. However, on a domain 
wall with de Sitter geometry, the trace anomaly acts like an extra contribution to the tension 
which permits a self-consistent de Sitter solution to the equations of motion. The Euclidean 
version of this is a spherical domain wall separating two balls of AdS. The radius R of the 
domain wall is given by p4j 

R 3 



P V I 2 8vr/ 3 
The metric in each bulk region is pure AdS: 
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ds 2 = l 2 (dy 2 + sinh 2 ydQ, 
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(7.4) 



(7.5) 



for < y < uq. The domain wall at y — yo, where yo is given by R = I sinh y^. 

The RS model can be interpreted using the AdS/CFT correspondence as four dimensional 
gravity coupled to a Yang-Mills theory with an ultraviolet cut-off |^6|, 24]. The Yang- Mills 
theory is two copies of the Af = 4 U(N RS ) super Yang-Mills theory with N RS given by 



G 5 
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(7.6) 



We shall refer to this dual Yang-Mills theory as the RS CFT in order to distinguish it from 
the theory on the domain wall. The Newton constant in four dimensions is given by the RS 
value G4 = G$/l. The four dimensional dual of the RS model with a U(N) CFT on the domain 
wall is four dimensional gravity coupled to both the RS CFT and the U(N) CFT. These two 
CFTs are rather different in that the former has an ultraviolet cut-off (so its effective action 



13 c 



See 1 4£ ] for an explanation of why this term is required. 
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does not behave as p 4 logp at large p) whereas the latter does not. The effective action of the 
RS CFT is proportional to iV| g , while the effective action of the other CFT is proportional to 
N 2 . This implies that the effects of the RS CFT should be negligible when iV ^> Nrs. This is 



confirmed by expanding equation |7.4| in powers of N/Nrs- At leading order, one recovers the 



four dimensional result 2.19. Note that N 3> Nr$ implies R^> I, i.e., the domain wall is large 



compared with the AdS length scale. 

7.3 Brane-world perturbations 

The RS model is a short distance modification of gravity. For length scales much greater than 
the AdS length I, four dimensional gravity is recovered. However, at shorter distances gravity 
becomes five dimensional. One might expect this to lead to a characteristic signal in the CMB. 
This turns out not to be true when the Yang-Mills theory is included on the domain wall. 
The reason is simple: at short distances, the matter contribution to the graviton propagator 
completely dominates the contribution from the four or five dimensional Einstein-Hilbert action. 
One might think that this effect is peculiar to our model of anomaly driven inflation, and would 
not occur in other models of inflation. However, any model has to take account of the Standard 
Model, which contains a large number of fields. These matter fields will suppress small scale 
metric perturbations in the same way as our Yang-Mills theory. 

We shall illustrate this effect explicitly by calculating the scalar and tensor graviton propa- 
gators for anomaly driven inflation in the RS model. Our method will be the same as above, i.e., 
we shall calculate the propagators in Euclidean signature and analytically continue to Lorentzian 
signature. The initial quantum state of perturbations is defined by the boundary condition of 
regularity on the Euclidean solution. In the RS case, this condition of regularity extends into 
the bulk. 

This work is an extension of our previous paper [2~3|, which contained the first rigorous 



derivation of cosmological perturbations in RS cosmology. However, in that paper we only 
discussed tensor perturbations and did not include the finite R 2 counterterm. Here, we shall 
include this counterterm and also consider scalar perturbations. Our method involves integrat- 
ing out metric perturbations in the fifth dimension. For alternative approaches to brane-world 
cosmological perturbations, see |T7|, p| [49], 1501 



The metric perturbation on the domain wall can be decomposed as in subsection |3.1| , giving 
a scalar ip(x) and a tensor 6ij(x). Correlation functions of these quantities can be calculated 
by integrating out the bulk metric perturbation, as explained in |f24}| . This is done by splitting 



the bulk metric perturbation 5g into a classical part 5go and a quantum part <5g'. The classical 
part is the solution of the linearized Einstein equation in the bulk that is regular throughout the 
bulk and matches onto the metric perturbation at the domain wall. The quantum part vanishes 
at the domain wall. Performing the path integral over 5g' gives some determinant Z that we 
shall not worry about. The classical part simply contributes the bulk action evaluated on shell: 

d[5g\e- Sbulk[5s] = Z e- Sbulk[5so] . (7.7) 

We conclude that the effective action governing metric perturbations on the domain wall is 

&eff — 2Sf, u jfc[5go] + Sbrane- (7.8) 
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The factor of 2 is necessary if we regard S^ik as the action of just one of the bulk regions. 
Sbrane is straightforward to compute using our result for W, equation |3.73| . The bulk metric 
perturbation 5go can be obtained from the results of section [3] by replacing I and G by I and 
G5. It follows that the bulk metric perturbation is transverse traceless, and the scalar ip arises 
from a perturbation in the position of the domain wall in Gaussian normal coordinates. Sbulk 
can be obtained from equations |3.69| , |3.7U| and |3.71| since the bulk action in the RS model is 
exactly the same as the bulk action for the AdS/CFT correspondence. 

From S e ff one can read off the metric propagators. The Euclidean scalar correlator can be 
written in position space as 
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The tensor correlator is 

(e ij (x)e i , f {x')) 
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F(p,y ,a,P) = e yo smhy Q 



'I'M 



+ 4coth?/o - 6 + *(p) 



fviVo) 

+ 2f3p(p+l)(p + 2)(p + 3) -Aap(p + 3). 



(7.12) 



Recall that tjq is defined by R = lsinhy . We have used equation [T.4| to write l 3 /G^ in terms 
of R. f p is defined in equation |3.56| . Equation [7.12| was derived in [23] but the term involving 
a was not included. In comparing our propagators in the RS model with those of the four 
dimensional model, we first render them dimensionless by dividing by R A . 

The scalar correlator for the RS model is very similar to that of the four dimensional model, 
as given by equation |3.83| . The only difference is the ^-dependence of the tachyon mass m 2 . 
As yo — > 00, the four dimensional value is recovered. This is to be expected since, in this limit, 
R/l — > 00, which implies N/Nrs ~^ 00 using equation fT4 . We have already discussed how the 
RS corrections are expected to be negligible when 3> Nrs- Note that as yo increases from 
to 00, m 2 just changes monotonically by a factor of 2/3. 

The analytic structure of the RS tensor propagator is very similar to the four dimensional 
case. There is always a pole at p = 0: this is the massless graviton of the RS modeQ- Other 
poles behave as discussed in subsection [0| 

The tensor propagator appears to exhibit more interesting dependence on y . The first term 
in equation [7.12| arises from the gravitational part of the action, so this is where differences 
between a RS model and the four dimensional model show up. As y$ — > 00, the first term tends 
to p 2 + 3p + 6, in agreement with the four dimensional result (equation 3.82|) . For very small yo, 



*This pole was mistakenly identified as gauge in [E4[ 
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Figure 4: F(p, y , 0, 0) plotted against p. The lower curve on the left (upper curve on the right) 
is for R ^> I, when four dimensional gravity is recovered. The other curve is for R = I, when 
the RS corrections might be expected to be large. However, they clearly have very little effect. 



the first term is p + 6. If yo is held fixed but large then the first term grows quadratically with p 
as p is increased but eventually becomes linear for sufficiently large p, corresponding to gravity 
becoming five dimensional at short distances. Thus the difference between a RS model and four 
dimensional gravity might be expected to show up in 1/p behaviour in the tensor propagator at 
large p, rather than the usual 1/p 2 behaviour. However, this neglects the effects of the matter 
fields, which are given by the other terms in equation |7.12| . At large p, ^(p) grows like p A \ogp 
and completely dominates the first term. Therefore, at large p the tensor propagator behaves like 
(p 4 logj9) _1 irrespective of whether one is considering a RS model or four dimensional gravity. 
The differences between the RS model and four dimensional gravity are drowned out by the 
damping effect of matter fields at short distances, rendering them unobservable. 

RS corrections are expected to be important at distances of order I. If we take R = I then 
all the tensor harmonics have wavelengths smaller than /, not just the large p ones. Therefore, 
one might expect RS corrections to be important at small p for such a small domain wall. 
Surprisingly, this turns out not to be the CclS6, clS shown in figure [|. This surprising behaviour 
can be understood in the four dimensional dual picture. Taking R = I corresponds to iV 2 pa 
6ANr S , so the matter on the domain wall still dominates the effect of the RS corrections. The 
RS corrections would be expected to be about as important as the matter on the wall when 
Nrs ~ N, which corresponds to R pa 0.46/. In other words, the RS corrections only become 
large when the entire domain wall is smaller than the AdS radius. 

One might worry that introducing a cut-off into the matter theory would spoil the damping 
at large p. However, if we did have a momentum cut-off A then we would need AR 3> 1 in order 
for field theory to be valid during inflation, as is always assumed. It therefore seems appropriate 
to take A ~ m p i, which corresponds to p ma x ~ A ^> 1. Figure ^ shows that the matter fields 
dominate the propagator even for quite small p, so introducing a cut-off would have little effect. 
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8 Conclusions 



There is now good observational evidence suggesting that the early universe underwent a period 
of inflationary expansion. Most theoretical models of inflation involve a scalar field rolling down 
its potential. The simplicity of such models is attractive but they have several serious problems. 
All these models require contrived initial conditions - no explanation is given of why the scalar 
field was initially displaced from the minimum of its potential^- Secondly, in order to obtain 
sufficient inflation and small CMB fluctuations, the CMB potential has to be highly fine-tuned. 
Finally, models of scalar field driven inflation usually disregard the effect of the large number 
of other fields in the universe. It is usually argued that the effect of such fields rapidly becomes 
negligible during inflation. However, as we have seen, this is not necessarily true because the 
trace anomaly of matter fields provides an additional contribution to the cosmological constant 
constant during inflation. 

In this paper, we have argued in favour of Starobinsky's model of trace anomaly driven 
inflation [15j as an alternative to scalar field driven inflation. In Starobinsky's model, the trace 
anomaly supports a de Sitter phase of expansion which is unstable, but can be long lived. This 
model is better motivated from the point of view of initial conditions since quantum cosmology 
predicts that the de Sitter universe can nucleate semi-classically via a four sphere instanton fi% . 
We have seen that this model admits a second instanton. This can probably be interpreted in a 
similar way to the Coleman-de Luccia [53| instanton, i.e., as describing the semi-classical decay 
of the de Sitter phase via nucleation of a pair of bubbles, each containing an open inflationary 
universe. Owing to the lack of an analytic solution for this instanton, we have concentrated on 
the four sphere instanton in this paper. 

During the de Sitter phase, particle masses would have been small compared with the space- 
time curvature so matter fields would have been classically conformally invariant. Moreover, we 
observe a large number of fields today and supersymmetry predicts that there should be many 
more, so the large N approximation is justified in studies of trace anomaly driven inflation. This 
leads to a very attractive way of calculating the effective action of matter fields during the de 
Sitter phase, viz the AdS/CFT correspondence. Using AdS/CFT, we have presented the first 
calculation of scalar and tensor metric propagators for trace anomaly driven inflation, taking 
full account of the back-reaction of matter fields. 

In order for the de Sitter phase to be unstable, it is necessary for the coefficient d = 
aN 2 / (1Q7T 2 ) of the V 2 i? term in the trace anomaly to be negative (in our conventions). We 
therefore included a R 2 counterterm in the action to control this coefficient. We also took 
account of the other curvature squared counterterms. We demonstrated that the amplitude 
of long wavelength metric perturbations could be brought within observational bounds at the 
expense of fine-tuning of iV 2 |a|. This fine-tuning is no worse than required in scalar field driven 



inflation, and agrees with the results of Vilenkin [17|. In fact, the amount of tuning required may 



be much less than for scalar field driven inflation. A more detailed treatment of the slow-roll 
phase would be required to verify this. 

One might worry that introducing a R 2 counterterm into the action would lead to observa- 
tional consequences for, say, solar system physics. However, the effect of this term in flat space 
is just to introduce a scalar field whose mass is of order m v ij {NyJ—ot). Even though \ct\ is very 



is 



Quantum cosmology can answer this question, but only for very contrived false-vacuum potentials EHl p4 
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large, this mass is still much too large to lead to observable effects today. For example, taking 
N = 10 and a of order 10 9 gives a mass of order 10 _6 m p ;. 

Our tensor propagator exhibits interesting analytic structure. We have shown that ghosts 
can be removed without fine-tuning, although this introduces a pair of complex conjugate poles. 
Such poles were studied long ago and found to correspond to violations of causality. We have 
seem that this causality violation occurs on a time scale Rj y/—a, where R is the Hubble time. 
This time scale is much smaller than R when |a| is large enough to bring the amplitude of 
metric perturbations within the observational bound^. 

At large p, the tensor propagator exhibits the behaviour first discovered for flat space by 
Tomboulis [22|, namely suppression of metric perturbations by matter fields. This suppression 
does not involve fine-tuning, as required for suppression of long-wavelength perturbations. The 
matter fields make the tensor propagator decay like (p 4 logp)~ 1 at large wavenumber p. This 
behaviour would be expected whenever the large N expansion is valid. Since we observe a large 
number of matter fields, we have argued that this suppression should occur even if inflation 
were not driven by a trace anomaly. This implies that matter fields damp out the effects of any 
short distance modifications of gravity (such as extra dimensions), rendering them unobservable. 
We illustrated this effect by comparing the propagators for trace anomaly driven inflation in 
four dimensions and in a Randall-Sundrum model. At large p, the tensor propagators are 
indistinguishable and at small p they only differ when the entire domain wall is smaller than 
the radius of curvature of the fifth dimension. 

There are many directions in which our work could be extended. For example, our use of 
AdS/CFT has restricted us to a strongly coupled theory. However, we have argued that our 2- 
point functions are independent of the Yang-Mills coupling. Dependence on the coupling would 
be expected to show up in higher order correlation functions of metric perturbations. This 
implies that these higher order correlation functions would not be determined by the 2-point 
functions, so the spectrum of CMB fluctuations would exhibit non-Gaussianity. 

In the Einstein static universe, the strongly coupled Yang-Mills theory exhibits a confine- 
ment/ de-confinement transition at a certain temperature, corresponding to two different bulk 



solutions in the AdS/CFT correspondence ||31|| . One might therefore wonder whether there is a 
bulk solution different from pure AdS which could have a spherical boundary with an 0(4) sym- 
metric metric. If so, then one might have a phase transition in a cosmological background. This 
does not appear possible. To see this, assume that the 0(4) isometry group of the boundary 
implies a corresponding 0(4) isometry group in the bulk (we are thinking of a cut-off CFT, cor- 
responding to a finite boundary). Birkhoff's theorem then implies that the bulk is (Euclidean) 
Schwarzschild-AdS. However, in order for the instanton to have spherical topology, the orbits 
of the 0(4) group have to degenerate at two points on the instanton (the poles) and this is not 
possible if the bulk is Schwarzschild-AdS except when the mass parameter vanishes. In other 
words, there is a unique solution (pure AdS) in the bulk and therefore no phase transition. This 
bulk solution corresponds to a deconfined phase of the Yang-Mills theory (this is evident from 
the overall N 2 factor in the Yang-Mills effective action). 

When one has a choice between several cosmological instantons, one usually argues that the 
instanton with the least Euclidean action is preferred, on the basis that this instanton would give 



16 Even if the time scale for causality violation were the Hubble time, it is not clear that this would contradict 
cosmological observations and such violations would certainly not be observable in the laboratory. 
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the dominant contribution to a gravitational path integral. Instantons which are saddle points, 
rather than local minima of the action, would not be viewed as satisfactory. These instantons 
would possess negative modes, corresponding to directions in field space along which the action 



decreases. Such instantons have been extensively discussed in ]5q| , where it was argued that 
they may be interpreted as describing quantum tunneling in an existing universe, rather than 
creation of a universe from nothing. Since we have found two instantons, it would be interesting 
to examine whether they have negative modes. This could give support to the idea that the 
double bubble instanton describes an instability of the de Sitter vacuum. 

Any discussion of negative modes presupposes the existence of a gravitational path integral. 
This is not well-defined even for Einstein gravity since it is well-known that the Euclidean 
gravitational action is unbounded below. In our case, the presence of the R 2 counterterm with 
a negative coefficient appears to make matters even worse. However, it is known that Einstein 
gravity coupled to a R 2 term can be rewritten as Einstein gravity coupled to a scalar field 
so the situation is probably no worse than usual. 

We have emphasized that there are two instantons in the Starobinsky model. However, there 
is also a third, namely flat space viewed as the infinite radius limit of the four sphere. This has 
infinitely negative Euclidean action. It might therefore be necessary to invoke the anthropic 
principle to explain why an inflationary universe is nucleated rather than an empty flat universe. 
The situation is analagous to false vacuum decay |53| , 54| , for which the instanton describing 



nucleation of a universe in the true vacuum state has lower action than the instanton describing 
nucleation of a universe in the false vacuum state. Clearly there is plenty of scope for future 
work on understanding the quantum cosmology of the Starobinsky model. 

Our approach was based on decomposing the metric perturbation into scalar, vector and 
tensor representations of 0(5), or 0(4,1). This made the AdS/CFT calculation relatively 
straightforward, but means that our results are only directly applicable to the initial de Sitter 
phase, although we have argued that the amplitude of metric perturbations should not increase 
during the slow roll phase. In order to produce a detailed fluctuation spectrum that could be 
compared with observation, it would be necessary to do a calculation based on a decomposi- 
tion into scalar, vector and tensor represenations of 0(4) (assuming a closed universe). If the 
AdS/CFT calculation could be extended to perturbations around a Euclidean background with 
a general 0(4) invariant metric then, by analytic continuation, one could calculate how the 
metric propagators evolve during the slow-roll phase. The perturbations spectrum at the end of 
inflation could then be used to predict the detailed spectrum of temperature fluctuations in the 
CMB. An 0(4) approach would also be necessary to investigate the double bubble instanton. 
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